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Abstract 

General implications of existence of non-perturbative scales and hadronic 
sub-structure for high energy processes are discussed. We propose that 
the dependence of the cross section of qq dipoles on their size d should 
deviate from when d becomes comparable to substructure scale. Then 
we discuss Kharzeev-Levin pomeron model [Q], based on ladder-type di- 
agrams with scalar resonances (scalar tttt or a and the scalar glueball 
Go). This channel is truly unique, because instanton- induced attrac- 
tive gg interaction leads to unusually small sizes and strong coupling 
constants of these states, supplemented by unusually large mass scale. 
Mo ~ 4 GeV, of the transition boundary to the perturbative regime. As 
pomeron is a small-size object by itself, these resonances may play a spe- 
cial role in its dynamics. Furthermore, we use more realistic description of 
the scalar gluonic spectral density without free parameters, and slightly 
modify the model to get correct chiral limit. We conclude that the non- 
perturbative part of the scalar contribution to the soft pomeron intercept 
is A = .05 ± 0.015, with comparable contributions from both a and Go. 

1. Significant progress of non-perturbative QCD has been mostly related 
with approaches based on its Euclidean formulations: lattice simulations, semi- 
classical theory based on instantons etc. During the last decade we have learned 
a lot about correlation functions and their spectral densities, hadronic wave 
functions and form-factors. Dramatically different features of different channels 
pointed out in ^ was confirmed and studied in details Q : below we discuss one 
of the most striking cases, the gluonic J^*^ = 0++ channel. 
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However, little of this progress has so far contributed toward understanding 
of high energy processes. True, it is difficult to translate many of those tools into 
Minkowski space. Elsewhere we will report some semi-classical calculations 
aiming to bridge this gap, while we start this Letter with more general discussion 
of some qualitative ideas related hadronic substructure to high energy scattering. 

The most important lesson we learned is that the non-perturbative objects 
in the QCD vacuum (and inside hadrons) are not some shapeless soft fields, 
with typical momenta of the order of Aqcd l/^^-^^j as it was assumed in 
70's. Instead we have semi-classical small-size instantons and very thin QCD 
strings. The instanton radius peaks around the p ~ l/3fm (|§l,for recent lattice 
data see and general review Q). String energy (action) is concentrated in a 
radius of .2 fm (.4 fm) in transverse directions Both are small compared to 
typical hadronic size, suggesting a substructure inside hadrons. 

A snapshot of parton distribution in a transverse plane inside the nuclcon 
should look like indicated in Fig.(|l]), for different x regions. These parton clus- 
ters originate from "scars" in the vacuum, being perturbed by external objects 
- valence quarks and strings, and therefore they must have the same transverse 
dimensions. One expects that these images of constituent quarks, diquarks and 
strings should be best seen at some intermediate x, before hadrons become 
black disks at very small x (high energies). The non-perturbatively produced 
sea quarks supposed to be more concentrated inside the constituent quarkaH 




while the string is supposed to be gluonic. A diquark cluster is also believed 
to be an instanton effect Finally, strong < qq > modification inside the 
nucleon should result in additional small density of sea quarks and gluons filling 
the whole disk (shown by light grey in Fig.l). 

If one prefers to use the language of hadrons rather than fields, existence of 
two distinct components can be viewed as being due to two different scales for 
glueball and pion clouds, respectively. However, using hadronic description in 
transverse plane (or t-channel) is probably not very useful, because such com- 
plicated and coherent field configurations as instantons and strings can hardly 
be discussed well in this way. 

These qualitative ideas were discussed in literature for long time. Con- 
stituent quarks as clusters were discussed e.g. in and "scars" of two strings is 
behind Nachtmann-Dosch model jlo| of high energy hadron-hadron scattering. 
But they are still mostly ignored by high energy practitioners, who think about 
partons as being randomly distributed inside the hadronic disk. 

How can we tell whether such substructure really exists experimentally? The 
simplest process we have is deep- inelastic scattering (DIS). In the target frame, 
it can be viewed as a scattering of a dipole-like qq objects with variable size 
d « 2/Q, where Q is the momentum transfer^. The small-d dipoles measure 

^This concentration should be enhanced for the polarized part of the sea. Sea quarks are 
found to be polarized opposite to valence quark (and the nucleon), as the instanton-based 
mechanism demands H. 

^This estimate works better for longitudinally polarized virtual photons, while for trans- 
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Figure 1: A snapshot of the parton density in a nucleon at different x, in an 
impact parameter plane. 

only the average gluonic field in the target^, but dipoles with d comparable to 
substructure scales p indicated above should show a nontrivial behavior of their 
cross section a{d) on d. pQCD predicts that a{d)/(P is constant at small d, but 
above some critical value d > dc ^ p we expect this ratio to drop, because such 
large dipoles start to miss the "black spots" in the target. 

Phenomenologically, when HERA data were translated into such dipole cross 
section ^ similar behavior a{d) was indeed found, and at the right scale 
for d. In principle however, two different effects can be responsible for it. One, 
emphasized in |p^ , [l^ especially for very small x is saturation of the cross 
section, when the unitarity bound (or "blackness") is reached. Our idea suggest 
similar behavior of a{d), but even at larger x ~ 10~^ where the spots are 
still "grey" and the dipole cross section is not large enough to need shadowing 
corrections. In fact, the particular parameterization used in have exactly 
this feature for all x. Further work is needed here to understand real magnitude 
of both effects, due to shadowing and substructure, respectively. 

Another source of information about parton correlations in transverse plane 
is diffraction. Clearly an inhomogeneous distribution we advocate enhances 
it, as compared to the homogeneous disk: the blackness of spots in higher, 
and there are more edges at which diffraction may take place. Studies of how 
diffraction depends on at HERA is therefore very important. One possible 
parameterization is hard-plus-soft pomerons p^ , which places the boundary 
between soft and hard contributions at ^ 4 — lOGeV^, oy d ^ dc we speak 
about. 

More generally, like it happened for form-factors at similar some time 
ago^, after new round of works we may be forced to reassess where exactly DIS 

verse ones large d tail is more significant, see [Q. 
''For exact definition see 

* Although dependence of form-factors roughly follow perturbative power counting 
rules, their absolute magnitude significantly exceeds the pQCD predictions. Non-perturbative 
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is truly perturbative. The previous paradigm - the leading twist dominance 
down to Q as low as ~ 0.5 GeV^ - seems now oversimplified. 

Hadron-hadron collisions are of course much more complicated. Studies of 
inelastic diffraction in 7rp,pp, pA and collision of two nuclei definitely show 
very large 0(1) fluctuations of the total nucleon cross section. It is clear that 
it is completely incompatible with the picture of a grey parton disk filled with 
multiple independent partons, with dozens of degrees of freedom involved. Only 
very few degrees of freedom may drive those fluctuations. One obvious suspect 
is the distance between constituent quarks, or the length of the string in Fig.l. 
Another is intermittent blackness of each constituent quark, presumably re- 
lated to "twinkling" character of the field strength distribution in the instanton 
vacuum. 

The pomeron itself is known to be a small-size object in transverse plane, 
as can be inferred from the fact that most of the t-dependence of pp scatter- 
ing is explained by nucleon form-factors, and also from a' ~ {2GeV)~^. It 
means when we see diffractive scattering of two nucleons, those are related with 
diffraction of their small partsj^ 

The history of pomeron goes back 40 years: it still works very well^, new 
discoveries are being made (such as existence of a pomeron polarization vector 
but its microscopic dynamics still lacks theoretical understanding. pQCD 
promised to explain the "hard pomeron" and although recent calculation 
[ pT| of the next-to- leading correction may put it in doubt, hopefully there will be 
a way out ^ Non-perturbative approaches include incarnations of the old multi- 
peripheral model (e.g. ||l^): but they neither provide clear cut explanations 
of why particular hadrons should be used, nor give convincing quantitative 
predictions. They also have no connection with Law-Nussinov 2-gluon exchange 
model, which seems to be a very natural starting point explaining constant (not 
growing with s) part of the cross section. 

New model for growing cross section and soft pomeron have been recently 
proposed by Kharzeev and Levin (KL) [Q. It includes (i) a ladder made of 
two t-channel longitudinal gluons, as in perturbative approach; (ii) while the 
s-channel "rungs" of the ladder being replaced by production of scalar physi- 
cal states. Their main motivation was to apply some known non-perturbative 
matrix elements, such as gluon-to- tttt transition near threshold. Furthermore, 

approaches, such as instanton-based calculation of the pion form-factor ]22) , are in quantitative 
agreement with data. 

^Except at extremely very high energies, when the whole disk becomes black. 

^Similar things happened before: let me just recall an optimistic example which may be 
unfamiliar to many high energy physicists. Free energy for high temperature QCD has negative 
0{as) correction of modest magnitude: but higher order ones show much larger corrections 
of different sign. However (unlike the pomeron intercept) one can calculate this free energy 
non-perturbatively on the lattice. The result is about 15% below free gas, close to the O(os) 
term. All high order corrections apparently canceled out! Furthermore, there are indications 
that some resummed or improved perturbation theory exists, in which those cancellations 
happen explicitly. 
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they inserted "realistic" scalar spectral density instead of perturbative one be- 
low some mass Mq, suggested a schematic model for it and estimated a result- 
ing value for the pomeron intercept]^ A w 0.08, close to the phenomenological 
value Q . Their input was however very limited and therefore they have treated 
gluonic scalar spectral density in a very schematic way. It explains correctly 
qualitative features of the result (e.g. its dependence on the number of colors 
and flavors N^Nf), but one may question the accuracy of the resulting num- 
bers. Below we (i) provide further motivation for the KL approach, and also 
(ii) improve on their crude schematic model in several respects. We explain 
why one should modify the expression (]l|) for the pomeron intercept, and get its 
meaningful chiral limit. We use other information about scalar spectral density, 
with different glueball parameters. We found that quark-related (pion) contri- 
bution to the pomeron intercept A is reduced, and the glueball one is enhanced 
compared to KL numbers, making them comparable at the end. 

2. To motivate the approach, let us start with the following question: How 
using hadrons in a ladder diagram can be consistent with the statements made 
above, that the pomeron exchanges take place between small parts of the col- 
liding hadrons, such as constituent quarks or strings? 

Well, it depends: different hadrons have different sizes! For example, mul- 
tiple papers use chiral Lagrangians and pions propagating inside the nucleon]^ 
It follows from chiral Lagrangian that it can be used for momenta up to the 
so called chiral scale p < « 1 GeV. Note: it is the small pion size which 
matters here, not its mass. We will have similar situation for glueballs below. 

With this in mind, we can move to the next question: Why is it reasonable to 
single out scalar 0++ gg channel? (Apart of the fact that we know few related 
coupling constants.) Because both its prominent resonances - the qq state a and 
the scalar glueball we call Go - are very small-size, a is the pion's brother and its 
interaction is covered by the same chiral scale. Remarkably, Go is even smaller, 
with Rgo ~ .2/m jlj, It is the smallest hadron we know, setting a record 
of its kind. It should be possible to construct some effective Gq Lagrangian, 
applicable below some momentum scale Mq. (For first attempts to build it see 
for discussion of the magnitude of Mo see below.) The qualitative reason 
why it is so compact is very strong instanton-induced attraction in the scalar 
channel due to small-size instantons. As shown in H, |^, in this case diluteness 
of the instanton vacuum (p/R)^ (1/3)* (where R « ^ fm is instanton mean 
separation) is compensated by classical enhancement factor {Sn'^ / g'^ (p))"^ ^ 10^. 

Other channels do not have this feature. For example, another gg channel 
one may think of is tensor 2++. However instanton field do not have such 

^ A enters the total hadronic cross section energy dependence as O'(s) ~ s^. 

* Well known ultimate model of that kind was suggested by Skyrme: in it a nucleon is 
made out of pions entirely. As noted by Witten, large Nc makes the nucleon static and pion 
field classical. Nevertheless, the ijjv^x ^ condition is still needed to justify the model, 
and because this parameter is not really large, Skyrme model cannot be very accurate at any 
Nc. 
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component, and so they do not act in it. Consequently the tensor glueball 
has normal hadronic size, R2++ ~ .8/m according to and it would be 
meaningless to consider its propagation inside a nucleon. As experience with 
quark vector channels (where the situation is similar) shows, in such cases we 
have multi-hadron spectral density dual to perturbative one down to low scale 

M2++ « Mq. 

KL have shown that the scalar channel contribute the following "non-perturbative 
part" to the pomeron intercept: 

A = ^ / ^ iPr>.ys{M^) - /-*(M^)) (1) 

where pphys{M'^), p'"'^'^{M^) are physical and perturbative (gg cut) spectral den- 
sities respectively. The integrand in (|l|) is non-zero only for M < Mq because 
at M > Mo pQCD works and two spectral densities become identical. 
3. Let us discuss properties of scalar spectral density Pphys 

(M^), first in glu- 

odynamics and then in QCD. 

We use the same normalization of the probing operator as KL 



so that its perturbative spectral density is 



1 b'g^ {N,'-1)M^ 

PpertiM) = — — (3) 



4096 

The spectral density should obey the low energy theorem 
dM2 



M2 



■[Pphys(M2) ~ ppert(Af^)] = -4 {OmO)\0) (4) 



As alsomphasized in , this relation has historically provided the first indication 
that there should be large non-perturbative scale in scalar gg channel. Later 
more direct and quantitative assessment of this effect was proposed [Q, based 
on small-size instantons. 

In gluodynamics p{M) is dominated by the contribution of the scalar glue- 
ball Gq: the lightest (and therefore stable) particle of this theory. Although 
phenomenologically its assignment to the observed scalar resonances is con- 
fused by mixing with qq resonances and is still under debate, both multiple 
lattice works and the instanton model point toward Mq^ — 1.5 — 1.7 GeV. 
As mentioned already, even more important is its small size, which leads to a 
remarkably large coupling constant to gg current^ which according to |Q is 

Ao =< 0\g^GlJGa >« 16. ± 2 GeV^ (5) 



^We remind the reader that the units in the gluodynamics is traditionally defined by setting 
the string tension to be the same as in QCD. 
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Substituting glueball contribution to the spectral density pca = (V (327r^))^AQ(5(M^- 
Mq^) one finds the foUowing contribution to the pomeron intercept 

Ago « -03 (6) 

This is not yet the complete non-perturbative part: one still has to subtract the 
"missing" perturbative contribution for M < Mq. 

As in we determine Mq from a "duality" sum rule[^ 

[PGo~ Ppert)^^ =0 (7) 

which ensures that correlator at small distances is not changed by changing fr'om 
Ppert to Pphys- Solviug it for Mq, ouc gct^ A/q ~ 4/, GeV. For a boundary 
between hadronic and partonic descriptions it is unexpectedly high scale indeed. 

The resulting contribution of "missing perturbative states" below Mq leads 
to negative contribution to A of about -.01. In total, we got ^giuodynamics ~ -02, 
about twice the value estimated by KL. 

4. Now we return to the real world with light quarks, and consider the sigma 
(or tttt) contribution. The major input of the KL paper is the tttt coupling at 
small M. It follows from the scale anomaly for chiral Lagrangian |Q 

P.. (8) 

which is larger than ppert because there is no [^j. The KL schematic model 
assumed that p{M) — p^^^^M) for all M < Mq (see the dashed line in Fig.(^)). 
But this assumption cannot be true for larger masses, because the pions are 
known to interact strongly in this channel, forming the famous scalar a reso- 
nance. 

Since we know its parameters, low energy tttt contribution into the corre- 
lation function in question can be easily reconstructed. Consider the isovector 
vector (p-meson) channel, for which the spectral density is well known from 
e -|- e— collisions and t lepton decay. In this case the tttt contribution at the 
threshold is trivial: the pion coupling is just the pion charge. Due to pion at- 
tractive interaction, the spectral density grows from threshold, till it reaches 
the peak - the p-meson. Its magnitude can therefore be fixed from the "vector 
dominance" , using the normalization to the M=0 point and known p meson 
width. 



^"Rather than from the low energy theorem (j^. With the simple sharp cutoff we use one 
cannot satisfy both, so we select duality because it closer to the integral we ultimately need. 



However the correlator calculated in ]2l[ is in exact agreement with this theorem. 

^^In KL paper significantly smaller number Mq 2.2 GeV was used for QCD. It is dual to 
only (T meson contribution, without a glueball. 
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Figure 2: Schematic representation of the spectral density, as i? = 
p{M) / ppert{M) versus the total mass M. The small circle at M=0 is the chiral 
Lagrangian prediction, the dashed line is the KL model, the solid line is our 
spectral density, which includes the contribution of the glueball, sigma meson 
and hadrons dual to perturbative gluons at large M. 



Similar "sigma-dominance" should work even better, because this resonance 
is very wide nia ~ To-. 

32 {AP - M^'^f + APrl ^ ' 

Naive integration from the threshold (2 m^) gives large contribution to pomeron 
intercept A]^"''"' w .09, where we have used M„ = .6 GeV, = .4 GeV. (Co- 
incidentally it is close to what KS got in their paper without sigma resonance, 
integrating till (their) Mq.) 

However, this large contribution cannot be correct, because it rest heavily 
on smallness of 771^^. In the chiral limit, when quarks and pions become massless, 
there is no threshold at 2m^, and the KL A simply diverges. 

To resolve this problem, we should know the distribution over the so called 
intrinsic parton transverse momentap]. If gluonic partons are indeed nothing 
else but expansion of classical field of an instanton, their transverse momenta 
are related to basic non-perturbative scale p, the mean size of QCD instantons. 
If so, one gets its right magnitude < >w (0.6 GeV)^ , and also predicts strong 
cut-offs, both at high and low momenta. The former is because the instanton 
is a finite-size object, leading to a form- factor ^ exp(ptp). The latter happens 
because the instanton field A'^ ^ V%^:^ has a vortex-like shape with changing 
sign, so that its projection to constant (or long-wavelength) field vanishesp]. 

^■^Deduced e.g. from pt distribution of Drell-Yan pairs, after perturbative effects due to 
extra gluon emission at large dilepton mass M are subtracted. 

In fact, there are experimental indications from diffractive dissociation cross section 
for transverse virtual photon is dominated by large pt, which probably implies that indeed 
xg{x,kt) ~ at small kt, see details in fccj]. 
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Therefore, gg collisions with small invariant mass are suppressed. 

Referring to quantitative calculation elsewhere ||2^, here we simply modify 
the KL expression of (|l|) with a logarithmic accuracy, introducing into the KL 
integral over transverse momentum of the exchanged gluons, J dk^/(kf + 
M^)^^ a finite cutoff fc™™. Now one of the factors l/M^ in (|^) is modified, and 

IStt^ /■ dAf2 (pp,^,(Af2)-pP-*(M2)) ^^^^ 



62 J M4 (M2 + (/fc™»)2) 

The unphysical divergence in the chiral limit is no longer there. Although 
sensitivity to the cutoff value is formally logarithmic, it still has significant effect 
on the intercept. With this modification, and (/c™*")2 for the pair to be equal 
to single parton < > mentioned, we find the a (or two-pion) contribution to 
be 

^a\realistic ~ -03 (1-^) 

It is now comparable to the glueball contribution discussed above. This outcome 
is in fact more natural than the KS numbers, 0.08 and 0.01 for a, Go effects, 
because they are is 0{Nj /N^) (for Nf — 2) and 0(1), respectively. 

5. Summary and discussion. Finally, combining all contributions together 
we get our final estimate of the non-perturbative contribution to pomeron in- 
tercept, resulting from hadronic ladder with M < Mo « 4 GeV to be 

A = .05 ±0.015 (12) 

now with our guessed uncertainties. 

How meaningful is this result? One still has to add to it the pQCD contribu- 
tion, which is ||l6|, |l7| with the region AI < AIq in the scalar channel subtracted. 
With current uncertainties, we do not know its value, and some readers may be 
disappointed at this point. However, the progress is not zero. First of all, the 
contribution we discuss in non-perturbative 0{a'^), enhanced by classical instan- 
ton effects 0(a7^) compared to perturbative result. Second, it is sufficiently 
close to the phenomenological value A — 0.08. Third, it can be experimentally 
tested by itself. 

Let us briefly indicate how it can be done. The discussed model claims that 
the ladders made of scalar resonances explains most of the growing part of the 
cross section, namely 6aNN « crjvAr(so) * 0.05 log (s/sq)- It means quantitative 
predictions about multiplicities of these scalar resonances a, Go, both in inelastic 
collision and their double diffractive production. Sigmas are wide and distorted 
by low pt cut we do not undertsnad well: so they are difficult to trace down. 
The glueball however shows up as relatively narrow resonance /o(1500), found 
in the double diffractive production. It has good experimental signature: large 
rjrj and decay modes. So its growing production with energy may even 
affect s-dependence of the rj/ir ratio. But the hottest thing to understand is the 
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azimuthal distribution of nucleons in double diffractive production: in fact the 
WA102 data (discussed in details in |26)) show surprisingly different distribution 
for "glueball-type" resonances (such as /o(1500) we discuss) from that for quark- 
based mesons. 
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